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ABSTRACT 

ON " 

Aims. By studying the photospheric abundances of 4 RVTauri stars in the LMC, we test whether the depletion pattern of refractory elements, 
seen in similar Galactic sources, is also common for extragalactic sources. Since this depletion process probably only occurs through 
interaction with a stable disc, we investigate the circumstellar environment of these sources. 
' Methods. A detailed photospheric abundance study was performed using high-resolution UVES optical spectra. To study the circumstellar 
^ environment we use photometric data to construct the spectral energy distributions of the stars, and determine the geometry of the circumstellar 
. - - i environment, whereas low-resolution Spitzer-IRS infrared spectra are used to trace its mineralogy. 

Results. Our results show that, also in the LMC, the photospheres of RV Tauri stars are commonly affected by the depletion process, although 
it can differ significantly in strength from source to source. From our detailed disc modelling and mineralogy study, we find that this process, 
as in the Galaxy, appears closely related to the presence of a stable Keplerian disc. The newly studied extragalactic objects have similar 
observational characteristics as Galactic post-AGB binaries surrounded by a dusty disc, and are therefore also believed to be part of a binary 
system. One source shows a very small infrared excess, atypical for a disc source, but still has evidence for depletion. We speculate this could 
point to the presence of a very evolved disc, similar to debris discs seen around young stellar objects. 
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1. Introduction 



With a distance of ~ 50kpc (Fea sjl999h . the Large Magellanic 
Cloud (LMC) is one of the best galaxies to study stellar evolu- 
tion. The known distance allows us to calculate the luminosity 
for stellar sources and it is the ideal laboratory to study physi- 



* This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. 



cal and chemical processes in an environment w ith a sub-solar 
metallicity of Z ~ 0.3- 0.5 Z dWesterlundl 19971) . Moreover, it 
is close enough to allow detailed studies of individual sources 
using large-aperture ground based telescopes. 

Here we focus on a sample RVTauri stars in the LMC. 
RV Tauri stars are pulsating evolved stars with a characteris- 
tic light curve showing alternating deep and shallow minima. 
They are located in the high luminosity end of the Population II 
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Cepheid instability strip dWallersteirj2002l) . The post-AGB sta- 
tus of RV Tauri stars was a long standing debate , but the d etec- 
tion of circumstellar dust around many objects (lJural 1986b . and 
the detection of extragalactic pulsators and their large derived 
absolute luminosities, were in line with the expected evolution- 
ary tracks of post-AGB stars. The first extragalactic RV Tauri 
stars in the LMC w ere discovered by the MACHO experiment 
dAlcocket al.lll998h. 



Revn iers & Van Winckell d2007l) and Revni ers et al.l (2007) 



performed a chemical study on two LMC RV Tauri stars, se- 
lected from those report ed by Alcock et al.l ( 1998 ). and found 
that, like in the Galaxy d Van Winckell 120031) . post-AGB stars 
are chemically much more diverse than previously anticipated. 
One of the stars, MACHO 47.2496.8, proved to be strongly en- 
hanced in s-process elements, in combination with a very high 
carbon abundance (C/0>2 and 12 C/ 13 C ^200). The metal- 
licity of [Fe/H] = -1.4 is surprisingly low for a field LMC 
star. The s-process enrichment is large: the light s-process 
elements of the Sr-peak are enhanced by 1.2 dex compared 
to iron ([ls/Fe] = +1.2), while for the heavy s-process (Ba- 
peak) elements, an even stronger enrichment is measured: 
[hs/Fe] = +2.1. Lead was not found to be strongly enhanced. 
The patterns can only be understood assuming a very lo w effi- 
ciency of the 13 C pocket (IBonacic Marinovic et al . 2007) which 
is created during the dredge-up phenomenon and the associated 
partial mixing of protons into the intershell. It was the first de- 
tailed study of the s-process of a post-AGB star in an external 
galaxy. 

Another object, MACHO 82.8405. 15, turned out to 
be chemically altered by t he depletion phenomenon 



( Rev niers & Van Winckell |2007). Depletion of refractory 



elements in the photosphere is a chemical process in which 
chemical elements with a high dust conde nsation temperature 



are systematically underabundant (e.g. iMaas et al .1 12005 



Giridhar et al.1 120051) . The special photospheric abundance 



patterns are the result of gas-dust separation in the circum- 
stellar environment, followed by re-accretion of only the 
gas, which is poor in refractory elements. The photospheres 
become deficient in refractories (as Fe, Ca and the s-process 
elements), while the non-refractories are not affected. The best 
abundance tracers of the depletion phenomenon are the Zn/Fe 
and S/Ti ratios because the elements involved in either ratio 
have a similar nucleosynthetic formation channel, but have 
very different condensation temperatures. Intrinsically Fe-poor 
objects have [Zn/Fe] and [S/Ti] close to solar, which is not the 
case for depleted objects. With [Fe/H] — -2.6, in combination 
with [Zn/Fe] = +2.3 and [S/Ti] = +2.5, in MACHO 82.8405. 15, 
there is no doubt that the depletion affected the photosphere 
of this LMC star very strongly. The very low abundances, as 
well as the clear correlation with condensation temperature, 
are shown in Figure HI which is a reproduction of the figure in 



Revniers & Van W inckel (2007). 



Photospheric deple tion is surprisingly common in G alactic 




12 3 
[3.6H8.0] 



Fig. 1. Colour-colour diagram indicating the Galactic post- 
AGB sources with discs (g rey circles) and e xtragalactic 



RV Tauri stars as presented by [Alcock et alj (1998) (black cir 



cles). The grey dots r epresent the LMC ob jects as found by the 
SAGE-LMC survey dMeixner et alj|2006l) . 



and very likely only fou nd around binary post-AGB stars (e.g. 
Van Winckel et al . 2009, and references therein). 

Also in their infrared spectra, the RV Tauri stars show 
unique spectral features. One of the best s tudied Galactic 
RV Tauri stars is AC Her dMolster et al] 1 19991) . Infrared stud- 
ies with ISO have shown very strong crystalline silicate bands 
from 10-50/im. Recent studies with Spitzer show that this 
strong crystallinity is commonly observ ed, also in other post- 
AGB binary sources dGielen et al.ll2008l) . 

Thanks to the efficient infrared detectors of the Spitzer 
satellite, very sensitive infrared observations allow us to 
probe for circumstellar dust, even around individual objects 
in external galaxies. The SAGE (Surveyi ng the Agents of a 
Galaxy's Evolution) Spitzer LMC survey ( Meixner et alj2 006) 
mapped the LMC, using all photomet ric bands of the Spitzer 
IRAC (InfraRed Array Camera, iFazio et aljl2004 and MIPS 



Rieke et al. 2004) 



(Multiband Imaging Photometer for Spitzer. 
instruments. This survey resulted in the detection of over 4 mil- 
lion sources. Thanks to the release of this database, we found 
that the LMC RV Tauri stars as discovered with the MACHO 
experiment in the visible, indeed have infrared excesses with 
very similar SED shapes as many Galactic post-AGB binaries 
(see Fig.Q]). Since only those LMC sources with strong enough 
optical fluxes were selected, some observational bias exists to 
stars where we see the disc more face-on, since edge-on disc 
would obscure the star too much. 

In this contribution we focus on the MACHO RV Tauri 
objects in the LMC as detected by the MACHO experiment 



post-AGB stars (e.g. Giridh ar et al.l 120051 : IMaas et al.1 12005L dAlcock et al.l 119981) . and we connect the chemical studies 



and references therein). In almost all depleted post-AGB ob- 
jects, there is o bservational evidence that a stable circumb inary 
disc is present dDe Ruvter et all2006HVan Winckelll2007l). The 
discs are very compact (e.g. Deroo et al. 20061: Derool 2007 ) 



based on high-resolution optical spectroscopy to the SED ener- 
getics and the infrared spectra obtained by Spitzer. Prior to this 
study, only 1 depleted post-AGB star in the LMC was known. 
Here we analyse the abundances of 4 more similar sources, test- 
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Table 1. List of stellar parameters for our sample sources. 



Name 


other name 


a (J2000) 


6 (J2000) 




l°gg 


[Fe/H] 


E(B - V) m 


Lm/L, 


L, 


P 






(h m s) 


(° ' ") 


(K) 


(cgs) 






(%) 


L 


(days) 


MACHO 79.5501.13 


J051418. 1-691235 


05 14 18.1 


-69 12 34.9 


5750 


0.5 


-2.0 


0.14 + 0.01 


60 + 3 


5000 ± 500 


48.5 




HV915 




















MACHO 81.8520.15 


J053254.5-693513 


05 32 54.5 


-69 35 13.2 


6250 


1.0 


-1.5 


0.27 ± 0.03 


0± 1 


4200 ± 500 


42.1 


MACHO 81.9728.14 


J054000.5-694214 


05 40 00.5 


-69 42 14.6 


5750 


1.5 


-1.0 


0.05 ± 0.02 


53 ±3 


4200 ± 500 


47.1 


MACHO 82.8405.15 


J053 150.9-691 146 


05 31 51.0 


-69 1 1 46.4 


6000 


0.5 


-2.5 


0.05 ± 0.01 


84 ±3 


4000 + 500 


46.5 



Note: The name, equatorial coordinates a and 6 (J2000), effective temperature r e fj , surface gravity log g and metallicity [Fe/H] of our sample 
stars. Also given are the total reddening E(B - V) tM , the energy ratio Lm/L,, the calculate d luminosity (comp uted by integrating the 
dereddened photosphere), assuming a distance of d = 50kpc, and the period P as given in dAlcocket all 19981) . 



Table 2. Log of the UVES observations and the obtained final 
S/N at a given spectral band. 



Star 


Exp. Time 


Wavelength 


S/N 




(sec.) 


(nm) 




MACHO 79.550 1.1 3 


3600 


375.8-498.3 


70 






670.5 - 1008.4 


100 


MACHO 81.8520.15 


3600 


478.0-680.8 


70 


MACHO 8 1.9728. 14 


3600 


478.0 - 680.8 


65 



Note: The S/N is measured in the middle of the spectral window 
covered. 



ing if depletion is also a common process in the RV Tauri stars 
of the LMC. As the distance to the LMC is known, we are able 
to discuss the evolutionary status of these extragalactic sources 
using their accurate position in the H-R diagram, which sofar 
has been impossible for Galactic post-AGB stars. 

This research was possible thanks to the SAGE-Spec inter- 
national program (http://sage.stsci.edu/). The goal of this large 
spectroscopic infrared program is to complement the wealth 
of photometric data from the SAGE photometric survey, with 
an extensive spectroscopic follow-up programme using the in- 
frared IRS spectrograph aboard of Spitzer (Kemper et al., 2009, 
submitted). The main goal of the survey is to determine the 
composition, origin and evolution, and observational charac- 
teristics of interstellar and circumstellar dust and its role in the 
LMC. A total of about 200 stars, in all stages of stellar evolu- 
tion, HII and diffuse regions were observed in Spitzer-IRS and 
MIPS-SED mode. 

2. Observations and data reduction 

2.1. Optical high-resolution programme 

The optical high-spectral -resolution data were obtained with 
the UVES dDekkeret al.ll2000b spectrograph mounted on the 
Nasmith focus of UT2 of the VLT, in February 2005 in visitor 
mode under mediocre sky conditions. 

To cover a wide spectral domain, two wavelength settings 
of the spectrograph were used: one with the dichroic, to use 
both arms of the spectrograph simultaneously (Dic2 860+437 
nm setting) and one centered on 580 nm using only the red arm. 



The detailed log of the observations is given in Table [2] The 
run was plagued with bad weather conditions, so for none of 
the three stars discussed here, we could obtain the full optical 
spectral window. 

The data were reduced using the dedicated UVES context 
in the MIDAS software package and include the standard re- 
duction steps of echelle data reduction. Spectral normalisation 
was performed by fitting polynomial functions through inter- 
actively determined continuum windows. For the lower S/N 
data, optimal extraction was used. Sample spectra are shown 
in Figures|2]and|3] 

2.2. IRS-spectroscopy 

In the follow-up programme SAGE-Spec to the SAGE LMC 
survey several post-AGB sources, showing RV Tauri-like be- 
haviour, were observed wit h the IRS instrument aboard the 
Spitz er Space Telescope dHouck et al.l l2004t IWerner etaL 
20041) . IRS was used in low-resolution SL and LL spectro- 
scopic staring modes. SL (A=5. 3-14.5 fim) and LL (^1= 14- 
37/im) spectra have a resolving power of R=A/ A A ~ 100. 
Exposure times were chosen to achieve a S/N ratio of around 
60 for the SL modes, LL modes have a S/N of 30. For a de- 
tailed description of the observations we refer to Kemper et al. 
(in prep.). Spitzer-IRS spectra for three of our sample stars are 
shown in Figure [5] 

The extraction was done from intermediate droop data 
products, version S17.0.4 up to S18.0.2, as delivered by 
the Spitzer Science Center. SMART reduction package tools 
dHigdon et al. 2004 ). and reduction tools developed by the 
FEPS (the Spitzer legacy program "Formation and evolution of 
planetary systems") team, were used for further extraction. For 
a detailed desc r iption of the different reduction steps we refer 
to iHines et al.l (12005b . The spectra are background-corrected 
and bad/hot-pixel-corrected and a fixed-width aperture is used 
for the extraction. The spectra are then combined and order 
matched. For the d efringing of the spectra the IRSFRINGE 
(Lahuis & Boogert 2003) package was used. Calibration spec- 
tral response functions were calculated, derived from standard 
stars and corresponding stellar models. 

In this contribution we limit the IRS spectral analyses to 
those objects for which we also have UVES data. 
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MACHO82.8405.15 
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Wavelength (nm) 

Fig. 2. Sample UVES spectra of the MACHO LMC stars. The spectra were normalised and offset for illustrative purposes. AC Her 
is a Galactic source which is a spectral analogue of the LMC sources. 
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Fig. 3. Sample UVES spectra of the MACHO LMC stars. The spectra were normalised and offset for illustrative purposes. AC Her 
is a Galactic source which is a spectral analogue of the LMC sources. The upper spectrum is from HR1017, a F5Iab star with 
solar composition. 
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Fig. 4. First and third row: The spectral energy distributions of our sample stars. Crosses represent photometric data, and the 
observed Spitzer spectrum is overplotted in gray. The solid line give the homogeneous disc model as discussed in Section [6] 
The dotted line gives the Kurucz stellar model. For MACHO 8 1.8520. 15 we only show the photometric data (crosses) and the 
appropriate Kurucz model (solid line). Second and fourth row: Below the corresponding SED we also show the abundance pattern 
[el/H] of the stars as a function of the condensation temperature. 
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Fig. 5. Spitzer-IRS spectra of MACHO 79.5501. 13, 
MACHO 81.9728.14 and MACHO 82.8405. 15. 



3. Photospheric abundance results 

One of the features of the spectra displayed in Figures[2]and[3]is 
the general weakness of metallic lines. The LMC RV Tauri stars 
are proven to be spectral analogues of the Galactic RV Tauri 
star AC Her, also shown in the figures. This Galactic RV Tauri 
star is considered to be one of the prototypes of its pulsation 
class, but it is also known to be a strongly de pleted object 
dGiridhar & Ferrolll989tlVan Winckel et al.lil998l) . 

Our analysis of the final produ c t spec tra is s imilar to those 
used in Van Winckel & Reyniers (2000) and iRevniers et al 



(20071). We performed a line id entification using the solar line 
lists of Th evenin dl989l Il990h . Weak lines, on the linear part 
of the curve-of-growth are the best tracers of the chemical con- 
ditions of the photospheres, as they are not saturated and are 
formed deeper in the photosphere, where possible non-LTE ef- 
fects are known to be weaker. For the chemical analysis, we 
only used lines with strength be low 150 mA an d accurate os- 
cillator strengths ( Van Winckel & Reyniers! 2000). 

As the stars seem chemically peculiar, we performed a 
relative abundance analysis using the Galactic RV Tauri star 
AC Her. The quantified abundance analysis was performed 
using the latest Kurucz model atmospheres of the appropri- 
ate metallicity (http://kurucz.harvard.edu/), and the LTE abun- 
dance analysis program MOOG (April 2002 version) written 
by Prof. C. Sneden (http://verdi.as.utexas.edu/moog.html). 

3.1. MACHO 79.5501. 13 
3.1.1. Photospheric model 

To determine the most appropriate photospheric model, we 
combined the classical spectroscopic methods with Balmer 
line-profile fitting and the constraints from the star being a 



member of the LMC. We refer to IRevniers et alj d2007l) for a 
more detailed description of the methods. 

The spectroscopic analysis assumes that the excitation of 
the Fe I lines follow a Boltzman distribution and by requiring 
that the abundances of lines from low and high excitation lev- 
els must be the same, a spectroscopic temperature is deduced. 
The gravity is determined by imposin g that the abundance of 



differ ent ions must be the same (e.g. lVan Winckel & Revniersl 
2000). For pulsating stars with large amplitude pulsations and 
associated shocks, this method turns out to be too sensitive to 
non-LTE effects. We therefore refined the effective temperature 
and gravity determination by fitting the Balmer 6 and y line 
profiles with synthetic models. As the object is known to be a 
member of the LMC, we assumed the distance to be 50kpc, 
and determined the absolute luminosity after dereddening of 
the colours. A surface gravity estimate is then deduced using a 
standard mass for a post-AGB star of 0.6 M . With this process, 
we find best values of = 5750 K, log g = 0.5 and a [Fe/H] 
abundance of —1.7. Typical errors are +250 K for the T e g, and 
0.5 for the logg. 



3.1.2. Chemical analysis 

The results of our quantitative analyses are given in Table [3] 
and displayed in Figure H] The spectral coverage of this ob- 
ject is not optimal as the intermediate spectral window of 
UVES with central wavelength of 580 nm was not observed. 
The results of our chemical study clearly confirm that the ob- 
ject is a spectral analogue to AC Her. Its Fe abundance is low 
([Fe/H] = -1 .8), even for the LMC, but this does not reflect the 
initial abundance. A clear correlation is found when depicting 
the abundance of a given element to its condensation tempera- 
ture (Fig. |4j: refractory elements are strongly depleted. Zn and 
S may be tracing the original metallicity which was likely to 
be around -0.6. This indicates that the object was originally 
slightly metal deficient compared to the overall metallicity of 
the LMC. 

Also in this LMC source, the depletion is found to be af- 
fecting the chemical content of the photosphere very strongly. 
For Sc, an element with one of the highest condensation 
temperatures, the abundance is about 250 times lower than 
what we would expect from a star of this initial metallicity 
([Fe/H] or igi na i = -0.6). Also Y and Ce, the only two s-process 
elements for which we found lines suitable for abundance de- 
termination, the underabundance is a factor 80 and 20 respec- 
tively. 

3.2. MACHO 81. 8520. 15 
3.2.1. Photospheric model 

The spectrum of MACHO81.8520.15 is of high quality which 
is reflected in the number of lines used in our chemical analyses 
and in a smaller line-to-line scatter. 

The star is found to be somewhat hotter than previous ob- 
jects. This is also seen directly in the spectra, more specifically 
in the ratio of line strengths between lines of ionised versus 
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neutral ions of the same species. Our best spectroscopic esti- 
mate is a r e ff = 6250 K with a gravity of log g — 1.0. 

3.2.2. Chemical analysis 

Also in MACHO 81.8520.15, the abundances are significantly 
lower than the mean metallicity of the LMC. Assuming the Zn 
abundance reflects the initial condition, a [Zn/H] = -1.4 makes 
MACHO 8 1 . 8520. 1 5 too a star of rather low initial metallicity. 

The abundance pattern shown in Figure|4]does show a gen- 
eral trend with condensation temperature. The main difference 
is that the C abundance turned out to be very high. The carbon 
abundance is determined on the basis of 8 lines which show 
little internal scatter, so the high C abundance is certainly sig- 
nificant. 

The low Zn abundance can be interpreted either as the ini- 
tial value, which means that the object increased its C abun- 
dance by a factor of 10, or it is the C abundance which reflects 
the initial conditions. In this case the depletion has been effi- 
cient to such low temperatures that even the Zn abundance was 
affected. 

3.3. MACHO 81. 9728. 14 

3.3.1. Photospheric model 

We obtained a good spectroscopic solution for the model at- 
mosphere parameters, showing that this star is a slightly cooler 
analogue of the previous object. The only Balmer lines covered 
are H^ and H a and these are affected by emission cores which 
make them unsuitable for model parameter determination. 

3.3.2. Chemical analysis 

The chemical pattern displayed by this star is depicted in 
Figure [4] This object has the poorest quality spectrum, and the 
line-to-line scatter in the abundances remains quite high. 

The abundance distribution is different in this object, since 
we find the Zn abundance to be very low ([Zn/H] = -1.22). 
Interpreting this as the initial chemical condition, this star is in- 
trinsically of low metallicity, even for the LMC. The abundance 
pattern is flat except for the refractory elements with the highest 
condensation temperature. These refractories show abundances 
which are significantly lower, down by up to -0.8 relative to the 
[Zn/H] abundance. Therefore, we interpret these abundances as 
affected by depletion. The effect has been marginal and only 
detectable for elements with the highest condensation temper- 
ature. 

The metallicity dependence of the depletio n phenomenon 



pletion patterns are abundant but only around stars where stable 
dusty discs ar e present (with some no ticeable exceptions like 
BD+39.4926 dVan Winckel et al.ll995b ). We the refore used the 



infrar ed catalogue released by the SAGE team dMeixner et al 



2006) to investigate the spectral energy distributions. We used 
the photospheric parameters from our UVES study and dered- 
dened the data until the match between the photometry and 
the Kurucz photospheric model was optimal. The SED figures 
are shown in Figure [4] It is clear that also around these ob- 
jects, which were selected on the basis of their light curve, cir- 
cumstellar dust is prese nt with specific colours typical of discs 



(De Ruvter et al. 2006). In the next sections we discuss our 
analysis of the spectrophotometric data we obtained for these 
objects. 



4.1. MACHO 81. 9728. 14 

For this star we unfortunately only possess spectral information 
in the 5 - 14yum region. The spectrum is unique in our sample, 
and is dominated by strong emission features at 6.2-7.6-8.6- 
1 1.3-12.8 and 14 urn, which can be identified as PAH emission 
(TJMenslbQOi and references therein). 

The C-C stretching and bending modes produce features 
with typical central wavelengths at 6.2 and 7.7 yum. The 8.6 fim 
feature is due to C-H in-plane bending modes and features 
long ward of lO/vm can be attributed to C-H out-of -plane bend- 
ing modes. 



Following the classification as described in iPeeters et al 



d2002h . which is based on the central wavelength of the main 
emission bands, the PAH emission bands are characterised as 
being of class A. This class A is usually linked to very pro- 
cessed PAHs residing in interstellar material, directly illumi- 
nated by a star. 

PAH bands are seen in the infrared spectra of several post- 
AGB sources, but surprisingly almost always together with fea- 
tures of silicate d ust species. Such mixed chemistry objects 



include H R4049 (Johnson et al. 1999 : iDominiket al 
Antoniucci et all 120051: Iffinkle et all l2007|v 



2003 



and HD 44179, 



the c e ntral star of the 'R ed Rectangle' ( Men'shchikov et all 



20021: ICohen et al.1 12004|) . In these stars the PAH emission 
is classified as of class B. Class B sources are mainly as- 



sociated with circumstellar material. HR4049 and HD44179 
are both binary post-AGB stars surrounded by an O-rich cir- 
cumstellar disc. The PAH carriers, however, do not reside in 
the disc of the se objects, but in a more recent C-rich bipo- 
lar outflow. In Giel en et al. I J2009h we describe some post- 
AGB sources where the less frequent class C PAHs are seen. 
Class A PAH features are seen in other post-AGB sources, 
has been found in the G alactic sample as well dGiridhar et al.l both post-AGB stars with evi dence for a carbon-rich chem- 
20051: iMaas et aL 2005), and may be confirmed here: at low istry, e.g. IRAS 16594-4 656 ( Garcia-Hernandez et al. 2006; 



initial metallicity, the effect is strongly reduced and often only 
visible for elements with the highest condensation temperature. 

4. Infrared spectroscopy 

Depletion was found to be very efficient in the RV Tauri stars 
for which we have UVES spectra available. In the Galaxy, de- 



van de Steene et al. 2008b as in mixed-chemistry sources, e.g . 
IRAS 16279-4757 dMolster et al.lll999l: uMatsuura et al.ll2004l) . 
In these two examples there is evidence for a dusty disc/torus 
and bipolar outflow in which the PAH carriers reside. 

If the PAH carriers reside in the disc then there has to 
be a formation process where carbon-rich species can be 
formed in an oxygen-rich environment. Such a scenario was 
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Table 3. Quantified abundance results for MACHO 79.5501. 13, MACHO 8 1.9728. 14 and MACHO 81.8520.15. 
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Note: The explanation of the columns is as follows: N gives the number of lines used; log e is the absolute abundance derived 
log 6 = log(N(X)/N(H))+ 12; £r M is the line-to-line sca tter; and [el/H] is the abundance relative to the Sun. For the references of the solar 
abundances nee ded to calculate the [el/H] values: see lRevniers & Van WinckeH fe007h . The dust condensation temperatures are taken from 
lLodderj J2003h . and references there in. They are computed using a solar abundance mix at a pressure of 10" 4 atm. The abundance table of 
MACHO 82.8405.15 can be found in lRevniers etail J2007h . 



proposed by lJura et al. (2006) to explain PAH features ob- 
served in the oxygen-rich giant HD 2335 17. The spectrum of 
MACHO 8 1.9728. 14 does not show any other features due to 
C-rich or O-rich dust species, such as SiC or silicates, although 
it cannot be ruled out that these features are present at longer 
wavelengths, for which no Spitzer-IRS spectrum is available 
unfortunately. Our chemical analysis of the photosphere shows 
that although the photosphere is affected by depletion, there is 
no evidence for carbon enhancement by a previous 3rd dredge- 
up episode. 



4.2. MACHO 79.5501.13 and MACHO82.8405. 15 




10 

Wavelength [yttm] 



MACHO 79.550 1 . 1 3 and MACHO 82. 8405 . 1 5 both show clear 
emission around 10 pm. This emission feature is mainly due 
to amorphous silicates, namely olivine (Mg2 X Fe2(i-. V )Si04) and 
pyroxene (Mg A Fei_ v Si03), which peak near 9.8 /im and 18pm. 
The observed shoulder near 11.3 pm is likely caused by crys- 
talline silicate dust, namely forsterite (Mg2Si04), the Mg-rich 
end member of the olivine family. We cannot exclude the pres- 
ence of other amorphous silicate species, but will focus here 
only on the olivine and pyroxene members. 



Fig. 6. Continuum-subtracted 
MACHO 81.9728.14 (black solid line) 
class A PAH emission features as described in iPeeters et al 
(2002). 



spectrum of 
Below we plot the 



MACHO 79.5501. 13 shows a strong slope in the contin- 
uum, and no strong features of amorphous or crystalline sil- 
icates at wavelengths longer than 10 pm are detected. Some 
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Fig. 7. The 10 - 14 - 16 - 19 - 23 - 33//m complexes of MACHO 79.5501. 13 (left), the Galactic mean spectrum (middle) and 
MACH O 82.8405. 15 (right ), continuum-subtracted and normalised. Overplotted in gray is the mean spectrum of the Galactic 
sources dGielen et al.ll2008l) . For the mean spectrum of the Galactic sources we also plot the comparison with mass absorption 
coefficients of forsterite (dotted), enstatite (dashed). In the 10 fim complex we plot amorphous olivine instead of enstatite, in 
dashed lines. 



emission might be seen around 15-19 and 23 fim. This shows 
that there is little cool silicate dust in the circumstellar environ- 
ment. The peak at 36 fim is probably a reduction residual at the 
end of the spectrum. 

In MACHO 82.8405. 15 some small features around 15 - 
19-23 and 33 /mi are seen. These features can also be attributed 
to the crystalline silicates forsterite and enstatite (MgSiOs). 
The 18 /mi amorphous silicate bump is detected, albeit not very 
strongly. 

4.3. Comparison with infrared spectra of the Galactic 
sample 

Both MACHO 79.5501. 13 and MACHO 82.8405.15 show a 
strong similarity to our sample of Galactic post-AGB sources, 
with similar spectral shapes and emission of amorphous and 
crystalline silicates. On average the Galactic sources do show 
a slightly higher fraction of crystallinisation. 

In our previou s study of the Galactic sample stars 
(Giele n et al.l 2008) we defined 6 different complexes (10 - 
14 - 16 - 19 - 23 - 33 fim complexes), and compared these 
to a calculated mean complex for the 21 Galactic sources. We 
will use this calculated mean Galactic spectrum to compare the 
observed features of these extragalactic sources. For this we 
subtracted a linear continuum from the different complexes and 
then normalised them. The comparison can be seen in Figure|7] 

The low signal-to-noise for these stars makes it hard to 
make a detailed comparison, but we can already see that sev- 
eral features appear in both Galactic and extragalactic sources. 
For both stars the 10 fim complex is very similar to the 
Galactic mean complex, although it seems slightly narrower in 
MACHO 82.8405. 15. The double-peaked structure due to the 
presence of both amorphous silicates and forsterite is apparent 
in both stars. 

Around 14/im we do find the enstatite 13.8/mi, and possi- 
bly the 14.4 fim, feature, but not the strong unidentified feature 
near 14.8 fim which is seen in the Galactic stars. 



For the Galactic sample we found the forsterite 16 fim to 
be shifted to shorter wavelengths. In MACHO 79.5501.13 we 
do not find strong evidence for this feature. The peak at 16 fim 
could be due to forsterite but note that the noise is large and 
enhanced by the continuum subtraction and normalisation. A 
strong feature at 15.5 fim is seen in MACHO 82.8405. 15. If this 
is due to forsterite the feature is shifted to even shorter wave- 
lengths in this star. In both stars there appear features around 
19 and 23 fim, but only MACHO 82.8405. 15 shows a strong re- 
semblance to the Galactic mean feature. The 23 fim complex in 
MACHO 79.5501.13 seems to be tilted to shorter wavelengths, 
instead of longer wavelengths. Whether this points to the pres- 
ence of a large enstatite fraction is unclear, since the other en- 
statite features seem to be only minimal. 

The 33.6/im forsterite feature is clearly present in 
MACHO 82.8405.15, albeit rather noisy, but can not be de- 
tected in MACHO79.5501. 13. 



5. Full spectral fitting 

For MACHO79.5501.13 and MACHO 82.8405 15 we per- 



form a full spectral fitting as already discussed in Gielen et al 
(2008). In short we assume the flux to originate from an opti- 



cally thin region, so we can make linear combinations of the 
absorption profiles to calculate the model spectrum. We allow 
for different dust species, grain sizes (from 0.1 to 4.0 fim) and 
dust approximations, such as Mie theory (Toon & Ackerman 
I l98l[) . Gaussian Random Fields (GRF, IShkuratov & Grvnko 
2005) and Distribution of Hollow Spheres (DHS, iMin et al. 
120051) approximations. We allow for two temperatures to de- 
scribe the spectrum. Unfortunately, the lack of strong emission 
features and the low signal-to-noise makes it hard to distinguish 
between different models, and the calculated^ 2 values are thus 
quite similar. Simil ar to what was seen in the Galactic sources 
(Gie len et al.ll2008l) . we find that the best fit is obtained using 
irregular grain compositions, and not a spherical description, 
such as applied in Mie theory. In Figure [8] and Table |4] we give 
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the result of our calculated best fit, consisting of small grains 
(0.1 - 2.0 /mi) in GRF dust approximation. 

For MACHO 79.5501.13 we find a good fit to the observed 
spectrum. The model continuum does not follow the strong 
downward trend, meaning that the modelled continuum tem- 
perature is too high. The 10 //m feature is nicely reproduced, 
as are the small bumps near 15, 19 and 23 /mi. The small fea- 
ture around 13.5 /mi is not explained by the model but could 
be due to enstatite, and was also seen in some Galactic sources. 
We find the model produces a significant 33.6 //m forsterite fea- 
ture, which is not observed and thus shows that the amount of 
cool dust is overestimated. 

A good fit is also obtained for MACHO 82.8405.15. Again 
the 10 /mi feature and features near 20 //m are very well fitted. 
There appears to be a rather strong feature at 15.5 /mi which 
is not reproduced in the model. Forsterite has an emission fea- 
ture at 16//m and for the Galactic sources we already found 
that this features se emed to be shifted to shorter wavelengths 
( Gielen et al.1 120081) . but the strength was generally well re- 
produced. The feature observed here however seems to be too 
strong to be explained by this shifted forsterite feature alone. 

6. 2D disc modelling 

To model the energetics of the obser ved RVTaur i stars , 
we use the disc model as described in Giele n et al. I d2007l) . 
nam ely the 2D radiative tran sfer passive disc model give n 
by lDuUemond et all d200lh and lDullemond & Dornimkl J2004 . 
The first part of this model consists of a 2D radiative transfer 
code which calculates the dust opacity and temperature. The 
gas component of the disc is assumed to be thermally coupled 
to the dust. The second part of the model is a ID hydrostatic 
equilibrium code which calculates iteratively the vertical den- 
sity profile. Stellar input parameters are the mass, luminosity 
and effective temperature of the star. For the total mass of the 
system we use a value of 1 M . We use the estimated lumi- 
nosity as given in Table Q] which corresponds to a distance 
d = 50kpc. The disc parameters consist of the disc size, to- 
tal mass, the different dust components and grain sizes and 
the surface-density distribution 2(r) = Y,o(r/ro) a , with a the 
power-law index and So a normalisation factor at reference dis- 
tance r<). The disc is made up from an astronomical silicate 
dust mixture with a gas-to-dust ratio of 100. We use a value 
-2.0 < a < -1 for the surface-density distribution, as expected 
in a disc environment. The inclination is kept fixed at 45°. 

The observed feature-to-continuum ratio can be used to de- 
termine the grain sizes in the disc, as small grains will produce 
much stronger feature-to-continuum ratios than l arger grains. 
For the Ga lactic sources we already modelled (IGielen et al. 
20071 l2009h . we opted to include an extra opacity source such 
as metallic iron and/or larger grains (up to 20 //m), but for these 
extragalactic stars we need to use considerably smaller grain 
sizes, ranging from 0.1 /mi to 5-10/mi. In the Galactic sources, 
we found grain growth to be very efficient in a disc environ- 
ment. This will produce dust particles with larger grain sizes, 
which will settle towards the midplane. This will cause the disc 
to be inhomogeneous, consisting of a disc made up from small 
grains and a cool midplane of larger grains. 



MACHO 79.5501.13 



0.06 



0.04 



0.02 



0.00 




MACHO 82.8405.15 




Fig.8. Best fits for MACHO 79.5501. 13 and 
MACHO 82.8405. 15. The observed spectrum (black curve) 
is plotted together with the best model fit (red curve) and the 
continuum (black solid line). Forsterite is plotted in dash-dot 
lines (green) and enstatite in dash-dot-dotted lines (blue). 
Small amorphous grains (0. 1 /mi) are plotted as dotted lines 
(magenta) and large amorphous grains (2.0//m) as dashed lines 
(magenta). 



For a few stars in the Galactic sample we possessed submil- 
limetre data showing a blackbody slope from 60 /mi to 850 /mi, 
pointing to the presence of extremely large grains, up to cen- 
timetre size or larger. The presence of such large grains in the 
disc can not be confirmed for these LMC stars since we lack 
data longwards of 40 /mi. 

The models are calculated using the canonical gas-to-dust 
ratio of 100. Since this va lue for the LMC may very well be 
higher (IGordon et al . 2003) because of the overall lower metal- 
licity, we tested how a higher ratio influences the disc structure. 
Increasing the gas-to-dust ratio means that dust formation and 
grain growth were less effective, and will thus give less pro- 
nounced dust signatures. We find that using a gas-to-dust ratio 
of 400, we need to increase the total disc mass by about a fac- 
tor 3 - 4 to get a similar dust mass, and thus to get a similar 
fit to observed infrared dust features. Changing the gas-to-dust 
ratio does not strongly influence the disc structure itself so the 
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Note: Listed are the^ 2 , dust and continuum temperatures and their relative fractions. The abundances of small (0.1 fim) and large (2.0 yum) 
grains of the various dust species are given as fractions of the total mass, excluding the dust responsible for the continuum emission. The last 
column gives the continuum flux contribution, listed as a percentage of the total integrated flux over the full wavelength rang e. The error s were 
obtain ed using a monte-carlo simulation based on 100 equivalent spectra. Details of the modelling method are explained in lGielen et al .1 
J2008h . 



Table 5. Results of our SED disc modelling. 
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Note: Given are the used grain-size distribution, the inner and outer 
radius (R m -R ont ), the total disc mass m for the homogeneous disc 
model and the surface-density distribution power law a. For 
MACHO 81.9728.14 we give both the results of the silicate disc 
model (S), and the carbon disc model (C). 



scaleheight of the disc is not very dependent on the gas-to-dust 
ratio. 



As noted in lGielen et al.l (120071 120091) the models are quite 
degenerate and equally fitting models with slightly different 
sizes, total masses and surface-density distributions can be 
found. Interferometric measurements, combined with the disc 
model, are invalua ble to constrain further the disc geometry 
dGielen et al . 2009). Unfortunately, the stars are too distant and 
too faint for current interferometric facilities, like the MIDI and 
AMBER instruments on the VLTI. 

6.1. MACHO 79.5501. 13 

For MACHO 79. 5501. 13 we find we need a very small outer 
radius (< 100 AU) to fit the observed SED. This could al- 
ready have been deduced from the spectrum, where the strong 
downward slope of the continuum pointed to the presence of 
little cool dust in the disc. A good fit was obtained by us- 
ing an outer radius of 30 AU and an inner radius of 10 AU. 
With these radii the inner rim has a temperature of 1200K, 
close to the canonical dust sublimation temperature for silicates 
(~ 1200 - 1500 K). The result of our modelling is shown in 
Figure |4] 



We find that at longer wavelengths the observed Spitzer 
spectrum falls more quickly than the model (even within the er- 
ror bars), which still overestimates the amount of cool dust. The 
observed infrared spectrum at the longest wavelengths is typ- 
ically less reliable so we would need longer wavelength pho- 
tometric observations to confirm the very steep slope at longer 
wavelengths detected in the infrared spectrum. 



6.2. MACHO 81. 9728. 14 

Modelling the SED of MACHO 81.9728.14 using the 2D disc 
model is less straightforward. The main problem lies in the fact 
that we only possess the PAH dominated Spitzer spectrum from 
5-15 /im, and no information at longer wavelengths. In other 
post-AGB sources the observed PAH carriers do not reside in 
the disc but in a more recent outflow, and so we cannot use the 
disc model to reproduce this part of the spectrum. Moreover, 
PAHs are small and subject to single photon heating effects, so 
they are not in equilibrium. 

If we do try to fit the SED with a standard silicate disc, we 
find that the model gives a strong 10 fim silicate feature, which 
should have been detected in the Spitzer spectrum. If the PAH 
features do come from the disc, then the silicate grains must 
be either very large, to avoid producing a strong silicate fea- 
ture, or very optically thick. Some small grains, however, have 
to be present to explain the observed near-IR excess. Another 
explanation would be that the disc is made up from another, 
probably carbon-rich, dust species. The chemical composition 
of the photosphere does not show evidence for a possible strong 
carbon enhancement. 

A disc consisting of 5 fj.m grains, and with inner and outer 
radii of 7 and 100 AU, proves a good fit to the observed spec- 
trum. A similar fit however can be found using a carbon- 
rich disc, with grain sizes of 0.1 - 10;um and a disc size of 
10-100 AU. The result of our modelling is shown in Figure|4] 
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6.3. MACHO 82.8405.1 5 

The 2D disc model also provides a good fit to the observed SED 
of MACHO 82.8405.15. For this star we find an inner radius 
around 17 AU. At this distance the inner rim reaches a tem- 
perature of 1000K, slightly below the dust sublimation tem- 
perature. The outer radius is less well constrained, but we find 
that it cannot be as small as in MACHO 79.5501. 13. In this star 
we see a clear 20 /mi silicate feature, which can only be repro- 
duced using an outer radius > 150 AU, with a value of 1.5 

for the surface-density distribution. The result of our modelling 
is shown in Figure [4] 

6.4. MACHO 81. 8520. 15 

The light curve o f MACHO 81.8520.15 was studied by 
Alcock et al. The classification of the light curve is un- 

certain but the period of 42.1 days makes it, together with the 
calculated stellar parameters, a possible RVTauri candidate. 
This star has a very low amplitude light curve with slight vari- 
ability in the de pth of its min i ma, b ut the period seems to be 
relatively stable. Alcock et al] ( 19981) suggest that it also might 
be classified as a long-period, low-amplitude W Vir pulsator. 

Looking at the spectral energy distribution of 
MACHO 8 1.8520. 15 (Fig. @}, we see that the star has a 
very small infrared excess, starting near 8 /mi. This is very 
different from the typical strong infrared excesses seen in 
RV Tauri stars. 

The small infrared excess shows that some circumstellar 
material is present, but it remains unclear what the exact loca- 
tion, geometry and origin of this material is. 

Interestingly the stellar photosphere is depleted. 
Photospheric chemical depletion is the result of a pro- 
cess which seems to occur only in objects surrounded by 
a dusty disc (see introduction), so we have at least indirect 
evidence that such a stable disc must be, or must have been, 
present around this object. 

The actual SED mimics that of a debris disc around 
younger stellar objects. These stars are surrounded by a very 
evolved circumstellar disc, which still contains trace amounts 
of dust and likely planetesimals of other macroscopic objects. 
These disc have been cleared of most of their original gas con- 
tent. Debris discs have b een found around main-sequence star s 
of all spectral types (e.g. lBrvden et al.ll2006HChen et all2006l) . 

As the photosphere is depleted, one could speculate the dust 
excess seen in MACHO 8 1 .8520. 15 as also coming from a very 
evolved and processed disc. Unfortunately, the lack of more 
photometric or spectroscopic data at long wavelengths makes 
it impossible to model the circumstellar environment of this 
star. 



7. Discussion 

The stellar parameters derived from high-resolution optical 
spectra, combined with ample photometric data, allow for a 
good estimate of the total reddening of the object. The high 
luminosities and the effective temperature prove the suspected 
post-AGB nature of these stars. All sample stars have effective 



temperatures ranging from 5750 to 6250 K and luminosities be- 
tween 4000 and 5000 L , assuming a typical LMC distance of 
50kpc. 

Our chemical study shows that, also in the LMC, the pho- 
tospheres of RV Tauri stars are commonly affected by the de- 
pletion process. The analysis of the Spitzer observations shows 
that, in the LMC as well as in the Galaxy, this depletion pro- 
cess is very closely related to the presence of a stable, dusty 
circumstellar disc in the system in which dust processing has 
been very active. The postulated long life-time of the dusty 
disc near the luminous star is a favourable circumstance for the 
gas-dust separation to occur, and for subsequent gas accretion 
dWaters et alll992h . 

For three sample stars we possess Spitzer low -resolution in- 
frared spectra. Two sample stars have spectra strongly resem- 
bling our previously studied Galactic post-AGB disc sources, 
with clear emission features due to amorphous and, in lesser de- 
gree, crystalline silicates. The atmospheric abundances, com- 
bined with the infrared spectral characteristics and the evidence 
for dust processing, show that RV Tauri stars in the LMC have, 
on average, very similar observational characteristics to the 
dusty RV Tauri stars in our own Galaxy. 

Based on our very limited sample here, the depletion pat- 
terns as observed in the LMC differ significantly from star to 
star, as was also observed in the Galactic sample. 

The most remarkable object in the infrared is 
MACHO 81.9728.14, which is dominated by emission 
peaks due to PAHs and does not show evidence for the 
presence of oxygen-rich dust species. In the photosphere, there 
is no evidence for an enhanced carbon abundance. The object 
shows some evidence of depletion, but only marginally so. 
The abundance pattern is flat with abundances around -1.0 
and -1.2 and with depletion affecting only the elements of the 
highest condensation temperature. With a [C/H] = -1.0 (based 
on two lines), the carbon follows this trend, and we interpret 
this total distribution as coming from an object with low initial 
metallicity. Despite this low initial metallicity, the dust excess 
is significant and the Lir/L* of 53% shows that the inner rim 
of the dusty disc covers a wide solid angle as seen from the 
star. The PAH carriers can be categorised as being of class A, 
which are the strongly processed interstellar PAHs. So far, it 
remains unclear whether the PAH carriers reside in the disc, or 
more likely, in an outflow from the central star. 

Also the abundance pattern of MACHO 81.8520.15 is re- 
markable. The star shows a very high [C/Zn] ratio and a very 
clear correlation of the underabundances with the condensa- 
tion temperature. If we interpret the high [C/Zn] value as due 
to depletion, the Zn abundance does not reflect the initial condi- 
tion, but must also be affected by gas-dust separation and sub- 
sequent accretion. This is remarkable as Zn is an element with 
one of the lower dust condensation temperatures. The alterna- 
tive interpretation, in which MACHO 81.8520.15 is an intrinsi- 
cally low metallicity object and carbon is enhanced by the 3rd 
dredge-up, is problematic: in such a scenario the s-process el- 
ements should be strongly enhanced as well, which is clearly 
not observed. 

Interestingly MACHO 8 1 .8520. 15 has a very small infrared 
excess, which seems to contradict the presence of a stable dusty 
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disc. However, as the star is depleted we postulate that a disc 
must have been present at some time during its evolution. The 
SED does show a resemblance to the debris discs seen around 
young stellar objects. In these debris discs, grain growth has 
already formed dust particles of considerable size (rocks and 
possibly even planetesimals), and the gas component has been 
removed. If there is a disc around MACHO 8 1.8520. 15 which 
is now strongly evolved, this would be the first time this late 
stage of disc evolution has been seen around a post-AGB star. 
The low Zn abundance would then indicate that the active gas- 
dust separation and accretion of the cleaned gas also takes place 
in colder regions of the evolving disc. 

MACHO79.5501.13 and MACHO 82.8405.15 are very 
similar in photospheric chemical composition, SED, and in in- 
frared spectral appearance to the strongly depleted Galactic ob- 
jects. The underabundance in some elements is a factor of 1000 
smaller than the solar value. The initial metallicity of both ob- 
jects is harder to recover, but the flat abundance distribution 
of C, Zn, S and Na suggests that both stars have an intrinsic 
metallicity of about -0.5, which is only slightly lower than the 
average value of the LMC. 

The spectral energy distributions of the sample stars are 
well modelled using a passive disc model. The inner radius 
is close to the dust sublimation radius for all stars, while the 
outer radius is poorly constrained, but a good fit is obtained 
with outer radii around 100-200 AU. We find, in comparison to 
most Galactic sources, that there is no need for an extra opacity 
source, such as metallic iron or larger grains, to reproduce the 
observed feature-to-continuum ratio of the dust emission. This 
could be an effect of the lower metallicity of the LMC sample 
which reduces the formation of metallic iron. Dust formation 
itself, and the subsequent processing, has been efficient enough 
to produce the large scale height needed for explaining the in- 
frared luminosity. Unfortunately, we lack submillimetre data to 
probe the component of large grains. Individual crystalline sil- 
icate bands do show some differences from the Galactic mean 
spectrum. To deduce whether this is a general trend in the 
LMC, and thus points to a different dust composition, grain size 
or dust model, we need to considerably enlarge the LMC sam- 
ple. This will be done in a future study. This extension of the 
sample will hopefully give a wider range in evolutionary phases 
in the post-AGB history of these sources, which we then hope 
to relate to other observed parameters, such as element abun- 
dances, disc parameters or infrared mineralogy. 

8. Conclusions 

Clearly the formation of stable dusty discs and the signifi- 
cant feedback of this disc on the central star is not exclu- 
sive to our Galaxy alone. Four out of five LMC RV Tauri ob- 
jects, for which we have high-resolution data, are found to be 
strongly affected by the depletion process in which the atmo- 
spheres became poor in refractory elements. Moreover, the in- 
frared colours and spectral data show that three sources are 
surrounded by a highly processed, stable disc in which dust 
processing has been efficient. In one source PAH particles are 
formed, while there is no other evidence for intrinsic carbon 
enhancement in the photosphere. This object is intrinsically the 



most metal poor object of the sample. For MACHO 8 1 .8520. 15 
we only found a small dust excess at 8 /an and we could inter- 
pret this as evidence for strong disc evolution. In this star, even 
the Zn abundance is affected by depletion which means that the 
gas-dust separation occurred at low temperatures. 

Overall the RV Tauri stars in the LMC display many char- 
acteristics of their Galactic peers. Whether these extragalac- 
tic stars also reside in a binary system remains unclear with 
the available data at hand. The low magnitude of these stars 
does not allow for a long-term radial velocity monitoring pro- 
gramme but, given the strong similarities to the Galactic bi- 
naries, a binary evolutionary channel is very likely needed to 
understand the RV Tauri stars in the LMC as well. 
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